%runpf(system7);
sysdata=system7;
busdata=sysdata.bus
branchdata=sysdata.branch;
branchi=branchdata(:,1);
branchj=branchdata(:,2) ;
busi=busdata(:,1);

branchz=branchdata(:,3)+branchdata(:,4)*1i;
tmp1=V(branchi)
tmp2=V(branchj)
Iij=(tmp1-tmp2)./branchz
Ia=abs(Iij)
Ib=angle(Iij)
% for i=1:length(Ia)

Iarel=Ia*100/110/sqrt(3)*1000